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ON FULLY DEVELOPED CHANNEL FLOWS: SOME SOLUTIONS 
AND LIMITATIONS, AND EFFECTS OF COMPRESSIBILITY, 
VARIABLE PROPERTIES, AND BODY FORCES 1 

By Stephen H. Maslkx 


SUMMARY 

An examination of the effects of comp feasibility , 
variable properties ) and body forces on fully developed 
laminar flows has indicated several limitations on 
such streams . 

In the absence oj a pressure gradient, but presence 
of a body force (e.g., gravity), an exact f ally developed, 
gas flow results. For a liquid this follows also for 
the case oj a constant streamwise pressure gradient. 
These motions are exact in the sense oj a Ornette 
flow. In the liquid case two solutions (not a new 
result ) can occur jor the same boundary conditions. 
An approximate analytic solution was jound which 
agrees closely with machine calculations. 

In the case oj approximately exact flows, it turns 
out that jor large temperature variations across the 
channel the efleds oj convection ( due to, say, a wall 
temperature gradient) and frictional heating must 
be negligible. In such a case the energy and momen- 
tum equations are separated, and the solutions are 
readily obtained. Ij the temperature variations are 
small, then both convection ejfects and frictional 
heating can consistently be considered. This case 
becomes the constant-property incompressible case 
(or quasi-incompressible case jor jree-convection 
flows) considered by many authors. 

Finally, there is a briej discussion oj cases wherein 
streamwise variations oj all quantities are allowed but 
only in such jorm that the independent variables are 
separable. For the case where the streamwise velocity 
varies inversely as the square root oj distance along 
the channel, a solution is given. 

INTRODUCTION 

Among all possible fluid flows, one of the most 
useful is the fully developed (i.e., independent of 

i Supersedes NACA Technical Note 4319 by Stephen 11. Maslen, 19.W. 


streamwise distance) channel flow. The flow is 
t akcn to be the motion generated by a const ant pres- 
sure gradient (the familiar Poiseuille flow) or by a 
body force (ref. 1). In either ease one usually 
considers only an incompressible or quasi- 
ineompressible flow' with fixed fluid propel ties. 
This is in marked contrast to the case of Couette 
flow where two parallel walls move with respect to 
each other. In such a case there is no need to 
limit oneself to a perfect gas or to any particulai 
variation of the transport properties (ref. 2 is a 
ease in point). 

The crucial difference between the Couette and 
the Poiseuille flows is that the former admits a 
stream wherein nothing depends on the stream- 
wise distance, while the latter requires that the 
pressure vary in the flow direction. Hence, to 
some small degree, atleast, the ot her fluid propert ies 
will also vary in that direct ion if t he state equation 
involves the pressure. The present study is ait 
examination of the general circumstances under 
which there can be a fully developed laminar flow 
past fixed boundaries. 

One special problem considered is an unusual 
situation found by Ostrach (refs. 1 and 3 to (i). 
He discusses the flow of a fluid in a two- 
dimensional channel under 1 1 ie influence of gi a\ it\ . 
Incompressible flow is assumed except as is ) e- 
quired to generate a varying body force, and the 
fluid transport properties arc assumed not to vary . 
Under these assumptions, the surprising result is 
found that there are two solutions to the flow in 
question for a certain range of values of the flow 
parameters. The first corresponds roughly t o the 
neglect of frictional heating, while the other is 
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Hour tin* (nontrivial) solution for homogeneous 
boundary conditions. 

There may be some doubt as to the stability of 
one of the solutions, presumably the second one. 
In any case, a question arises about the effect of 
considering a real fluid having variable viscosity 
and thermal conductivity as well as being truly 
compressible. 

Accordingly, the present paper treats the con- 
sequences of such generalizations. However, to 
reiterate, one serious restriction is made on all 
the flows considered herein: The flow is always 
fully developed, with the result that the effects of 
conditions near either end of the channel are 
ignored. 

ANALYSIS 

Consider a two-dimensional flow of a viscous 
compressible fluid acting under the influence of 
an axial body force such as gravity. Variable 
viscosity and thermal conductivity are admitted. 
The configuration is shown in the following sketch: 



1 he equations of motion are 

{p1 l )x-\-(pv) r=0 (1) 

p(uu A .+e« K ) +P X =- pf+2(nu x ) x +[ li (u r +v x )] r 

2 

— ^K«Y+»r)]x (2) 
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P (%+ <Wy) +I\ = i 2 2(nVy) y + [tl(Uy + V X )) X 

2 

^[m( w .yTiv)]j' (A) 

PC v {uT x -\-vTy) +P{U x +Vy) = (kT x ) x + ( kTy ) y 

+ p[2u x -\-Uy-\-v i x -\-2v$r-\-2u Y v x — (4) 

(Symbols aie defined in appendix A.) Consider, 
in addition, two possible forms of a state equation, 
one applying for a gas, and the other for a liquid: 

P=pliT (Gas) (5a) 

p = p[* —P(T—T)\ (Liquid) (5b) 

111 equation (5b), p is the (small) volumetric 

expansion coefficient, and p and T are reference 
values. The significant difference between the 
two state equations is that the second is independ- 
put of pressure. 

The boundary conditions on the channel walls 
are 

wX,0) =u(X,d ) =»(AT,0) =v(X,d) =0 

T(X,Q) = T Wo (X) 

T(X,d) = T at (X) 

where d is the distance between the channel walls. 
The temperature boundary conditions could, of 
course, be replaced entirely or in part by a heat- 
tiansfer condition but, for the purpose of this 
report, such a change is unimportant. 

Equations (1) to (ti) are sufficient to define the 
fully developed flow in a channel provided the 
viscosity and conductivity variations with tem- 
perature are known, and provided further that the 
forced-flow pressure gradient, if any, is specified. 
In seeking solutions of these equations for flows 
m a very long channel (i.e., fully developed), 
three approaches are considered: first, exact 

solutions entirely independent of distance along 
the channel (A); second, solutions approximately 
independent of A'; and third, solutions wherein 
the variables are separable. In each case the 
results ca i be expected to differ according to which 
of the stare equations applies. In t his connection, 
it is important to observe that the viscosity and 
thermal conductivity vary differently in liquids 
and gases. In particular, the viscosity rises with 
temperature for a gas and falls for a liquid. 
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The phrase “exact solution” should perhaps be 
defined. In this report it is understood to mean 
a solution which satisfies equations (1) to (t>) 
rigorously. However, no consideration is given to 
conditions near the ends of the channel. There 
are two relatively simple circumstances under 
which such exact solutions can be found. For a 
gas nothing can vary with X f not even the pressure. 
For a liquid this restriction is moderated to the, 
extent that only the gradient of the pressure need 
be independent of A r . This relaxed condition 
occurs because of the pressure-independent state 
equation (eq. (5b)) for a liquid. 

EXACT FLOW OF CAS 

If nothing depends on A", equation (1) to (5a) 
and (fi) become 

v=r() (7a) 

( y)y—pI (7b) 


Then equations (7) become 



P— Constant r 


Mr~* 


(ii) 


If the viscosity varies linearly with temperature 
(n,= l), then these equations are separated. If 
ri 2 is also unity, the solution is 


u= 



1) 


( 12 ) 



a (aJPPfV / V- 1\ 

t>\A RmIJ V 12 7 


03) 


P r = 0 (7c) 

(kTy^Y— — M u r (7d) 

P= P RT (7c) 

u( 0 )=u(<l)= 0 " 

T(0) = T W(i l (8) 

T(d) = T Wl ^ 


Assume that the viscosity and thermal con- 
ductivity vary as powers of temperature, and also, 
for convenience, change, the independent variable 

(Y). 

Thus, suppose 


and let 


or 


n=aT n i \ 
k=bT n * / 


i= — 1- 


2 r dY i 

Jo m/m « 0 

*. f-f- 

Jo M/Mv„ 


2 Y_ p/Jth; 

P J — 1 Mw 0 


(9) 


( 10 ) 


where the constant B is as yet undetermined. 


Finally, B can be found from the second form ol 
equation (10). Thus, 


2dT„ 2a/ a B^V 

IT = J , rd v=V"'+ r *J + \YbV4Rrt ) 


1 f /r/ '"'Vh 

240 Vd7’„J L\ 




-frm 


/here 


T 

1 m 


T + 7 T 

* JC 0 1 1 


= 1 (14) 

(14a) 


The quantity in brackets is essentially the 
parameter K defined in reference 1 This is al- 
ways positive. Under these circumstances equa- 
tion (Id) has only one real root, that root being 
such that 


0< 


BT m 

<1 K 


<1 


For example, for air under standard conditions, 
if d—^ and/ is gravitational acceleration, equation 
(14) yields BTJd T v = 0.904. A convenient 
standard for comparison of the present solution 
with more approximate results is given by the 
mass flow. This is 



BP r + 1 

2BT.J-1 


u&i}= 


/BT m y/pif(P\ 
\d T W J V 12 n m ) 


(15) 
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The corresponding result for incompressible flow 
with constant fluid properties, or for compressible 
flow with variable properties, but neglecting 
frictional heating, is 




For the case cited prior to equation (15), where 
BT m k\ 7^ = 0.904, tbe actual mass flow is about 
25 percent below the incompressible value. Un- 
fortunately, the mean flow velocity is about 800 
feet per second. This high velocity would proba- 
bly preclude the possibility of laminar flow even 
existing. If an example leading to slower flow 
were considered, the difference between the two 
results would have been small. This demonstrates 
that the variation of viscosity is unimportant. 

A final comment : If the viscosity is not assumed 
to vary linearly with the temperature, the mo- 
mentum and energy equations cannot be separated. 
As this circumstance (lack of separation) led 
Ostrach (rels. 1 and 4 to (>) to find two solutions 
for the flow rather than one, it is perhaps worth 
examining further. Define 




Then equations (11) become 


(16) 


_ a B 2 f Y / « // (// o — /i i + 1 ) \ ~ 0 ~ " i } ^ ( 1 - w i + w 2 } 
dAViioV b ) 


/>— Constant 


(17b) 


This completes the solution for the exact fully 
developed flow of a gas. There are two generaliza- 
tions which can readily be made. These involve 
the addition of a body force transverse to the chan- 
nel and the addition of heat sources in the fluid. 
The solutions are given in appendix B. 

EXACT FLOW OF LIQUID 

Here it is assumed that there is a pressure 
gradient such that, at least, P x and P Y are inde- 
pendent of A. No other A-dependence is ad- 
mitted. Then the system is that given in equations 
(7) and (8), except for t he .Y-momen turn equation, 
which can conveniently be written as 

(MHy)y T-T*)j 

I (Uv \ pf)-Ji3f(T-T*) ( 18 ) 

where 

T*= 1>x +p(H-J 7V ( j 9 

pPS 

It is worth observing here that, if the reference 
point is changed in equation (5b) (i.e., a new T), 
this has no effect whatever on the value of r T*. 
This is because the state equation (eq. (5b)) is 
really of the form p=A—BT , where A and B 
are fixed. In that case equation (19) is really 

rp* P X A f 

Bf 

lienee 2'* is a function only of P X) /, and the mate- 
rial. T len if rj and //, defined by equations (9) 
and (16 , respectively, replace y and T , the mo- 
mentum eq nation becomes 


For gases, one expects that /n<l and /> 2 >0; thus 
the right side of equation (17a) is a decreasing 
function of If. In this cast 1 there can be at most 
one solution of the problems. The argument goes 
this way. Suppose one solution is known. If a 
second solution has larger If } then by equation 
(17a), v vv is smaller. Hence, for a reasonable 
profile, m, is reduced. Then, by equation (17b), 
ff vr} is of lesser magnitude. Hence for a reason- 
able (‘use, If must also be small, which is a con- 
tradiction. 

It should be emphasized, however, that this case 
of no pressure gradient whatever has no connection 
with the work reported in references 1 and 4 to 6. 


"”'-(W)" <7 ’- r ' ) ,2o) 

while the appropriate energy equation is, again, 
equatioi (17b). 

Now the forcing term (the right side of eq. (20)) 
should be considered. According to its definition 
(eq. (16) ), If is an increasing funct ion of T. Hence, 
at least ‘or constant viscosity, the forcing term is 
an increasing function of If. If p and k are con- 
stant, then the forcing function is linear in H and, 
for this case, it has been shown (refs. I and 6) 
that two solutions occur under certain conditions. 
In the pi esent situation things are not that simple, 
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and the results depend on how ^ and k vary. 
For many liquids the conductivity varies only 
moderately over a fairly wide range of tempera- 
ture, while the viscosity may change severalfold. 
Two cases, water and liquid sodium, are illustrated 
in tables I and 11. In the present discussion the 


TABLE I.— VISCOSITY" AND CONDUC- 
TIVITY OF LIQUID SODIUM 



Conductivity, k, 

Viscosity 

. a, cenli- 

Turn per- 

watts/ (cm) (°C) 

poises 

at urc, 





T, °C 

Expert- 

1117 

Expert- 

Eq. (21) 


ment 
(ref. 9) 

IFto+t 

ment 
(ref. 9) 

(a) 

201) 

0. 815 

0.815 

0. 4.5ft 

0. 4.50 

300 

. 759 

. 700 

. 345 

. 3.50 

400 

. 712 

. 712 

. 284 

.286 

,500 

. 008 

. 008 

. 243 

. 242 

000 

. 027 

. 030 

.210 

. 210 

700 

. 590 

. 598 | 

. 180 

. 185 

800 

.547 

. 567 

. 165 

. 166 

900 




. 150 

. 150 

‘7’„= 150°( 

# = 158 (centipoises)(° 

0). 



TABLE II. — VISCOSITY OF WATER 
UNDER ATMOSPHERIC CONDITIONS 



Visc<» 

ity, ft, 


centi poises ; 

Temper- 
ature. r. 


- — 

°C 

Experi- 
ment 
(ref. 9) 

Eq. (20 
(a) 

0 

1. 79 

1.79 

20 

1.01 

.90 

40 

. 66 

.60 

60 

.48 

.45 

SO 

. 36 

.36 

100 

.28 

.30 


a T a = 20° C; 

/r= 36 (cent i poises) (°C). 


variation of conductivity is neglected. The vis- 
cosity can be written to good approximation as 


where s and T a are constants (for water «=0.36 
(centipoisq) (°C), and T a = 20° C if T is in °C). 
This expression is compared in tables I and II with 
experimental values. It is to be expected that 
(7N~ T a )>0 in the range where the fluid remains a 
liquid. Under these circumstances, equation (16) 
yields 

H=£-T{T+2T a ) (22) 


and equation (20) becomes 


/ B 2 pl3_fs\ 

lkTl+2sH(T*) 

\ 4<7l 

\T-71- 2sH{T) _ 


From equation (17b), it is seen that, if i/.j| ^ 
large, II must vary more or less parabolically 
upward across the channel. Then it follows (eq. 
(23)) that, when H—H(T*) has a large magnitude, 
the forcing term increases only slowly with //, 
while for small //— //(7 1 *), the forcing term is 
linear in II. The latter reduces the problem to the 
usual free-convection situation (ref. 1), while the 
former (large //— Il(T* )) approaches the usual 
Poiseuille case, wherein the forcing term is (‘ou- 
st ant. This circumstance at least restricts the 
range of flow parameters for which t wo solutions, 
as found in reference 1, can exist. 

Stated more explicitly, if (T— T^) / (T^r T a ) 
remains small, the viscosity is essentially constant 
and the system becomes that solved in reference 1 . 

Had the variation of conductivity been allowed 
for, a small modification of the foregoing argument 
might occur. If the conductivity drops as the 
temperature rises, the forcing term would move 
toward a more linear variation with H . 

In general, the solution of the system given by 
equations (23), (17b), and (8) is not simple. 
However, after two limiting cases arc discussed, tin* 
general case can be described. The first such cast 1 
is that of small frictional heating; the second is for 
small temperature variations, and therefore vis- 
cosity can be considered constant . 

SMALL FRICTIONAL HEATING 

First consider the case of small frictional heating. 
The formulation involving rj and II (eqs. (23) and 
(17b)) is not convenient for this case. Hence, 
consider equations (18), (7d), arid (8). The vis- 
cosity is defined by equation (21), and the con- 
ductivity is assumed constant. When the fric- 
tional heating is small so that the right side of 
equation (7d) is negligible, one obtains 

T=7\+j(T w - T w J (24) 

Then equation (18) can be integrated to give 

+ (12«i a.) (//) — 2«2«3 ()/)l (25) 
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wIltMV 


« 1 = 


T*-T* 


I k'q + I a 

*** T T 

1 W 1 1 fr o 


(26) 


« 3 = 


1 2a l a 2 jylfaHjgi ) +3 
(l+2o a ) 


and 7 * is defined by equation (19). From these, 
the wall sli ear is 


(m« r) r - o pj'&<l ( 7 ^ — T Wq ) (a3/12) 

( M M >' ) V ■■= tl~ pj ( / w | 7 ir 0 ) ( 0 : 3 / 12 Q!] 1 / 2 ) 

The net mass flow is 


(27) 


figure . The fluid is liquid sodium, and the 
temperatures of the two walls (100° and 900° C) 
differ enough that the viscosity varies by a factor 
of about 4. In spite of this there is no significant 
difference between the results for the two cases. 
That is, the effect, of variable viscosity is un- 
important even though the temperature variations 
are large. 

It is interesting to observe the case when 
T*—(T w +T Wi )l2, the average fluid temperature. 
Then a, = — 1 2, — 1 , and the shear is not only 

the same at each wall, but is the same in the 
constar t- and variable-viscosity cases. However, 
the velocity profiles differ slightly; and, while the 
mass flow is zero for the constant -proper tv case, it 
is not for variable properties. 

SMALL TEMPERATURE VARIATION 



rf )3 

360# 1 ^ Wq) 


^ I” } 1 )(1 +2q|)(l ^g) 

L l+2a 2 J 


These results can he compared with those for 
constant fluid properties. (This case is given on 
j). 10 of ref. 1.) In such a case, the temperature is 
again given by equation (24). The velocity 
distribution, wall shear, and mass flow are given, 
respectively, by 


■u= — 


p0f(l(T w 


12* 


7 "' J - a +-’«,) 


x [Ki“ 1 )(i +1+sh ')] (29) 


(nUy)r.o-=pf»KT r -T ir p [(l+:ia,)/6] "'l 

. r (:50) 

(pHy) y=a~ l T n ^) [ — (2 + 3a 1 )/6] J 

J’" (7„. -7’, n )’[(1+2a,)(l+2a 2 )] 

(31) 


In these last equations the viscosity has been taken 
as that corresponding to the average temperature; 
that is, 


2s 

F 'w x ~\~T v 0 -f-2 T a 


(32) 


A comparison of the results obtained for the 
cases of constant and variable properties is given in 


In the cast* whore the frictional heating is 
considered but where the temperature variations 
are small, the fluid properties can be considered 
const ai t. This ease lias been solved by Ostrach in 
some (I d ail (.refs. 1 and 6), by machine methods. 
However, another method of getting the same 
results is now presented that has the advantage 
of giving the parametric- dependence simply. The 
same method is applied later to the general case 
(large temperature variations). However, the 
justification of the procedure is most convincingly 
display mI by comparison with the aforementioned 
machine solutions. 


In this case of small temperature differences, it 
is again convenient to work from equations (18), 
(7d), and (8). If the temperature changes are 
small enough that ( r f—T W( )f(T Wi -\-T a ) is small 
everywhere, the viscosity and conductivity can be 
considered constant. 


The following dimensionless quantities are 
defined : 




(T r ,, VPPd' _ K(T-T*) 

T ~ { 1 } 


y os) 


9 Y 1 
v=z 7 -l 


where 


K= 


^ mMm 


CM) 
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is the parameter defined in reference 1. In terms f '(±\) — i) 

of the new variables, equations (18), (7d), and 

(8) heroin e r ( — 1 ) — Kj 1 () 



0 .2 A .6 .8 

Y/d 

Fkjukk 1. -Effect of property variation on velocity profiler. l\ t l()() c C; 7 T „ V 900° C; 7 T *, 0° C; ]i«]uici eodiuni 


513555-60 2 
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An approximate solution of this system can be 
found by iteration. The velocity profiles are 
usually parabolic*. Thus, suppose 

r (0 wi(i~r) (:*8) 

where A is an undetermined parameter. Then 
equation (3(1) yields, subject- to the boundary 
conditions, 

t=— [(m — 1)H _7? ('*9) 

If this is put into equation (35), the resulting- 
velocity distribution is 

'■•"--IDT 


This process could be continued, assuming con- 
vergence, but is stopped at this point. If equa- 
tions (39) and (40) are put into equation (30) 
and the result is integrated across the channel, a 
quadratic equation for the unknown parameter 
A~ follows. This is 


1 - / ct :i2r> r i— 1 ,r ri 

\ 80 S L 210 J / 


, 1,1 Of) 

+ 1 08,424 


X [A* ( //i I 1 

Heal solutions exist if 
210 

"/ 707 f 1 

\ 2,040 L (wil)-J 

' + \4>!b 40 [ 


1 — ^ 

4 ( m -1 1 )~_ 


— -0 < 4 1 ) 


</v( m + 1 ) 


‘210 


m 


( m T 1 ) 




(42) 


The boundaries defined by equation (42) are 
plotted in figure 2 for X> 0. Some limiting valued 
found in reference 0 by machine methods are 
shown for comparison. Agreement is excellent. 

The two solutions can be examined generally 
in the following manner. If /£(/« + 1) is moderate 
(say, in the range 0 to 40), then equation (41) 
yields, approximately, the two results 

A 2 ** [0.016K(m + l)] 2 


} 


(41a) 



Kiel ice 2.- Regions of existence of solutions of equations 
(35) l;< (37). 


Tlie nass flow and heat transfer to the walls 
are, res met ivoly, 


or 



k M Kim\ 1) 

tf<i 24 


Hi 


■5 


Y 

I 


Mk m I y =k m |( 7 r) r-ir ( I r.)r--rf] 

jmkj r.- T, 

:\K4 

2rM-J r T Wn -- T*) 

:\K<I 40 


(43) 




* r ( 44 ) 


The lird (smaller) solution is one in which there 
is negligible frictional heating and, hence, heat 
transfo \ while the second is quite the opposite. 
The se< oiul (*ase is, as is pointed out in reference 4, 
one of regenerative heating. There is a large* 
amount of heat transfer to the walls (eq. (44)); 
this heat is supplied by frictional heating of the 
fluid occasioned by large mass flow (eq. (48)) 



ON FULLY DEVELOPED CHANNEL FLOWS: SOME SOLUTIONS AND LIMITATIONS 


9 


.•uhI the resultant high shear. Notice that even 
when T w = 7 T » 0 = 7’*, so that the problem is an 
homogeneous one (K= 0 in rq. (37), although 
KKT*- T*) in oq. (33) is not zero), this second 
solution does not vanish. For that matter, the 
second solution is virtually independent- of K and 
w, provided K and m take on moderate values. 

T- he range of validity of these results is limited 
by the condition that (7 T — 7 , « !o )/(7 T » u + T*) be small 
so that the viscosity variation is negligible. In 
the case of the second solution (large A 2 ), the 
maximum of T occurs near the center of the 
channel, and thus equations (3S), (44), and (4ha) 
yield 

T m(U -T tH 200 7\-T* 
f ^ 7\ r ^ -f~ T„ 


If, for example, K^U), then 7^ — 7’* is limited 
to a few degrees. This implies also that the wall 
temperatures must be virtually equal. On the other 
hand, for applications of the first solution (small 
/l 2 ), the only restriction is that (7 \— r T W{ )!(T^ 
+ 7«) be small. However, if the viscosity is that 
for some suitably defined average, the error due 
to larger temperature variations should he unim- 
portant. This conjecture is based on the example 
discussed after equation (32). 

Before giving a numerical example, it is worth- 
while to examine the order of magnitude of the 
numbers one obtains in physical problems for the 
present case of small temperature differences. 
From equations (33), (34), (40), and (41 a), one 
lias, very roughly, 


u ( yj = 0 ) : 


[ k ( 7 w 

\ k u 


... fk(T~ — f*) 
- 2 ,) "••• 

\ v-K 





T(v 

O) 

For 

examp] 

It*, for 

water 

at 

7V 

-T* = 0. 

.25° C 

, and rn 

= 1, 

t he 

various 1 

rases are 





\ ( 

Small 



u(t> = 0), 

ft/sec 

2 r» 



Tin- 0) 

-r v »°c 

0 


This case corresponds to a channel width of 14.(> 
inches. A more general idea of the orders of the 
mini hers involved in the second solution can he 
obtained as follows. From equation (33), 

u 2 k ( 2 

T— T*~ n t ,40) 

However, from equations (39), (40), and (41a) 
the maxim ums of l and r for the second solution 
usually occur near 77 — 0 and have the respective 
values of about (>1 rt-nd 14. Hence, for the second 
solution, 


( t '/•*), 


- 7’*) J 

- ?’*) 
,-7-) 


(Small solution) 
(Large solution) 

(Small solu t ion) 
(Large solution) 


(45) 


Actually there is no combination of K and m 
such that equations (41), (39), and (40) yield 
f4» aJ c/T wai < 1. For water, 3k/n is of the order 
10 4 (ft/sec) 2 /°C, and hence if T— T* is 1° (\ the 
velocity maximum is 100 feet per second. For 
liquid sodium, 3k is about 10 7 (ft/sec) 2 /°C; thus 
a temperatures differeiiee of only 1° C corresponds 
to a maximum velocity of 3,000 feet per second. 
It therefore appears that, if the velocity is to he 
kept moderate to maintain laminar flow, the 
temperature variation must be small. Hence, the 
assumption of constant fluid properties is a 
good one. 

To compare the present method of calculating 
the velocity and temperature profiles with the 
more exact solutions obtained by machine methods 
(ref. 1), a single example is shown in figure 3. 
Remember that this is a constant-fluid-property 
situation. The example is that of water flowing 
in a channel 14.7 inches wide and for which 
T Wo —T*—A/10° C. The wall temperatures are 
20.0° and 20.1° C. This leads to A'=10. The 
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l’TcntK 3. Velocity and temperature profiles, m , 2: /;, 

10- T„ 7’* 1 10° ('; water. 

agreement is within 10 permit for the second 
solution and, not surprisingly, virtually exact for 
t he smaller one. 

Before ending the discussion ol exactly hilly 
developed Hows, it should be observed that the 
iterative procedure used to get solutions here can 
he applied in the other cases considered in refer- 
ences l, 4, and G, namely those involving wall 
temperature gradients and heat sources in the 
lluid. 

As is stated earlier, the present iterative proce- 
dure can be applied directly to the original problem 
wherein the frictional heating is considered and 
large temperature changes an* contemplated. I he 
solution for such a (low is given in appendix C for 
the ease of equal wall temperatures. The only 
difference from the ease just discussed is that 
some of the integrals are rather involved and the 
equation for the amplitude is more complicated. 


However, as is observed previously, if the velocities 
are to be kept moderate in the second solution, 
the temperature variation will he negligible. 

APPROXIMATE SOLUTIONS 

The results thus far presented have the beauty 
of being exact within tin* limitations of fully 
developed flow. However, several cases arise in 
which such a limited view is unacceptable. The 
simples* such east* is the How of a gas with a pres- 
sure gradient, the ordinary Poiseuille flow. An- 
other ease of some interest is that involving a wall 
temperature gradient. The extent to which such 
Hows can he considered fully developed is examined 
in what follows. 

Solutions of equations (1) to (G) are sought. 
Again, a long channel is assumed and end effects 
are neglected. In such a case, any gradients in 
the A"- direction (in the flow direction) must he 
small. Hence, write the variables as follows: 

r-T>{\ \ tu's) i ...| ^ 

p- p[p ( )(//H <5pi<T,7/) i . . d 
7 - T[T {) {y)-\ bT x {j\y) f . . .) 

p p(pn(//! -f <5pi (./',//) I • . d 

k — k [ hu ( // ) -(■ t ,// ) H - -d y <17) 

d7[ 7/,///) - 5 w-i (./*,?/ H • * d 

r= ( j-u\{) \ ; . . .] 

.i - Xj L 

u y/<i 

where l;L is small, d being the wall spacing and 1. 
being i length of flow, as yet undefined. The 
other two parameters, e and <5, are small hut un- 
related at this time. For a gas 5 must he at least 
as large as e for the state equation (eq. (5)) to make 
sense, while for a liquid the number 5 is determined 
b v the temperature boundary conditions. The 
barred quantities are, except for u, given param- 
eters chosen so that p 0f A*„, p (h and so forth, are of 
unit order. The value of u is initially unknown 
because there is no characterisin' velocity for in- 
ternal Hows of this kind. 

Had the term g(x) in the pressure (the first of 
eqs. (47)) been considered as a function of y also, 
then added terms would be introduced because of 
the //-momentum equation (eq. (H)). However, 
these are, analogously to the usual houndary- 
layer analysis, of higher order than what is re- 
tained in the other equations of motion. 



ON FULLY DEVELOPED CHANNEL FLOWS: SOME SOLUTIONS AND LIMITATIONS 


1 1 


If equations (47) are put into equations (1), (2), 
and (4) and only the dominant terms of each kind 
are retained, there follows 


Po( u i P r+i'i,,i/) 4 'PoV ? i+Pi,xMo ~0 ( 48 ) 


[pof%Wl.z+^l W 0.|/)] + ^ J ^+p/(po+5pi) 


} 


JJ2 Km0 W 0 ,y)y] (49) 


6pc r uT 


5 Pn 


- - [po ( v o 7 1 , j 4- r 1 7 o, y ) ] 4“ ~~j~ ( u \ . t + ?’i , y) 

Jf JjHP 

~ [(* 0 / 0 .y)y-\-S(ko 7 | , ?/ + A’j 7 (J y )J + (goM<j.j/) 


(50) 


Finally tlu‘ state equal ions (e(js. (5)) vied 


P=pllT 

Po/()= 1 

Po 1 1 "h Pi /o= (c/5 )(/(./’ ) J 

Po -l=— 0 7474-1) 

Pi — — / 377 i j 


> (Gas) 


74 74-1) j 

f (Liquid) 

artirj T I 


(51) 


(52) 


Now observe that the terms in the braces in 
equation (49) describe the driving forces of the 
(low and must therefore be of the same order of 
magnitude as the viscous terms (which contain 
the highest derivatives). The following four eases 
will bo considered in’ turn: 


Actually, it is only necessary that Pe/L<pf in order 
to have the body force matter. However, it can 
be assumed that equation (58) holds, and e can 
be discarded later if it is small enough. 

Then, using equations (51) and (52), the driving 
term in (‘({nation (49) is 


Y UiU') 4 p/(p<i)=p/ |j 4 — 1 J 


(Gas) 


(Liquid) 


(54) 


For a mixed flow, g z (x) — 0(1) for a gas or y T (.r)~ 
— l+0(/37 T ) for a liquid. A pure fr ce-conveetion 
flow might arbitrarily be defined as one for which 
<7*0r) = — 1, but, for convenient* e, any flow involv- 
ing body forces is henceforth referred to as a f roe- 
convection flow. 

Definitions of the barred reference values in 
equations (47) are all straightforward, except that 
for u, and a selection of values can readily be made 
a priori. However, u must be chosen such that u {) 
is of unit order, and there is no way of knowing 
ahead of time bow big the flow will be. Hence, for 
the moment, let us bog the question and define 
simply 

(',~=P.fFP 



■ Hotly 

Lurtre t em per- 

Case 

force 

attire variation 
( T n variable) 

, 

Yes 

Yes 

n 

No 

Yes 

m 

Yes 

No 

L\ P 

No 

No 


CASE I— To VARIABLE, FREE CONVECTION 

For ease I, terms of order 5 can be neglected as 
compared to corresponding zero-order terms. 
Thus, pi can be set equal to zero in equation (49), 
as can and T i u in equation (50). Assuming p n 
and (/Ax) to be of unit order, the body force* and 
pressure gradients must be of similar size so that 


where ('> is unknown and, for a gas, /JTs I. It is 
shown later that (\ is a number in tin* range 10 to 
50. Then equations (48), (49), and (50) can be 
written 

Po ( w i , T + 1' l . y ) 4“ P n , v l' \ 4" P J , j «0 = 0 f 5 6) 


(P( ) ?/()_ y 


viqV'-H+aA)]} l 


— [po('Wo'W 1 .r4 _ ^t'Wo.y)] (Gas) 

( Po W' 0 , y .) y + ( \ { 7 T Q I [1 4“ Hr (/ ) ] 7 1 } 

8u 2 ( '-> r . t 

= fj jj, lpo(‘MoM,, jr -rr l t£o t v)J (Liquid) 


r (57) 
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< h i\ a=4t |po(Mo7V 1+4,1 Tn “ )] +fy r ] 

x£po(Ho7'i.i+f>] 7 0, v) Mo~ (J x (j'j 

(fro7;.,),=^ ^3 ^[po(V7’i.x+Pi7o.»)l-r^- 
Xp 7 (Mo7’,,r+i’i 7’o jjy — 


J (Gas) 


(Liquid) 


The system whose solution is sought is given by 
equations (56), (57), (58), and (51) or (52) plus the 
boundary conditions (eq. (6)). Assuming that 
<j x (x) is given, the unknown dependent variables 
remaining are seven in number, w 0 , U\, V\ } Po? Pn 
7 T 0 , and 7\. With five equations and seven de- 
pendent variables, some further restriction must 
be made. The difficulty arises mainly in the terms 

and Pi. For a gas it appears p h u u and V\ are of 
the same order as 1\ and hence that the inertia 
term in equation (57) is of the same size as the 
convection term in equation (58). Hence, it 
follows that, in order for the solution to be de- 
termined, both of these terms must be negligible. 
The same result follows for liquids, although per- 
haps not so obviously. In this case (eq_ (51)), 
p, is very small, of the order of where is the 
volumetric expansion coefficient. However, for 
liquids t he viscosity is a very strong function of 
temperature, and thus p 1 = 0(7 7 i). lienee, there 
is no particular reason to assume that w i and 
are not the same size as 7\. Accordingly, if these 
terms matter, the present formulation is useless. 

For t hese terms to be negligible, two courses are 
open. One is to have everything ^-independent 
as in the exact solutions described earlier. The 
second is that the parameter (\hu*/JL$T he small. 
This is not simple. For example, if equation (55) 
is used and a channel 1 inch wide and having a 
characteristic length L of 10 feet is assumed, then 
under standard conditions and gravitational ac- 
celeration, 

( \&u 2 ijLirr=m,QQQ a/r 2 

-0,400 5/( 2 


(For water) ^ 
(For air) I 


( 50 ) 


Consider (\— 1. For the air case, 8 must be as 
large iu «, which is about 1/3,000 here. For water 
8 can be chosen by an applied wall temperature 
gradient. In either case it seems difficult to make 
C>8u 2 lfLfiT small, in fact, unless ( y 2 is a large 
number. 

If this question is ignored for the moment, it 
can he observed for a gas that, if 8^e (and 5 must 
be as large as e), the coefficients of the inertia or 
convection terms and of the dissipation terms are 


and 


r,bu l jfLa’l '= (' 2 tu 2 /fL=( 2 pu 2 /l J "] 


fiu~ in v 

FT ' o 


(pn 2 /P) 


y 


( 00 ) 


(where Pr is a Prandtl number), which arc both 
essentially the squares of a Mach number and are 
the same size if f r 2 = 1 . Hence, for the gas case, 
at least, the frictional heating must also be negli- 
gible because, as is shown later, f 2 >1. 

For x liquid it can be seen, by trying some rases, 
that t ie frictional heating must again be very 
small provided 8/e is not virtually zero. 

Now to return to the question of a value for ( 2 . 
Suppose that equation (55) is used together with 
its subsequent consequences. Then, for the sake 
of an example, let <j x (x) — — 1 and assume a gas. 
The differential equations become 


(jLtQ Wo, y)y — C 2 ( 1 ./ o) / 1 0 


(Mo.,), — 0 


The boundary conditions art* 

?/o(0)-i/o(l)-0 

7o(0) T v o 

/()( I ) — V , fl 

To -nuke the point about orders of magnitude, 
consider that mo and ko are proportional to the 
temperature (i.c., po— A* () = Yo). Then the equa- 
tions can be solved very read ily. Maximum speeds 
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have been computed for several cases and are the energy equation. The result in parametric 
shown in the following table: form is 


| HO | 

7V 0 7’ u |( 2 

i mar 

U 1 0 . 12 

!■» 3/2 .01 

1 2 | 03 

Then [-u-oUax is of unit order if C 2 is, say, 40. Thus 
for the conditions cited in connection with equa- 
tion (59), 

( i&iF/JLpT— 1 ,9505 (For water) 

or 

= 1005 (For air) 

which, particularly for the air case, can readily ho 
made small. In a case (such as eqs. (35) to (37)) 
where two solutions occur, the same conclusion 
about (\ follows, although the argument is rather 
tortuous. 

Finally, then, the inertia terms must he negligi- 
ble for the flow to be fully developed. In that 
cast' the thermal convection and frictional dissipa- 
tion are negligible. The equations for fully devel- 
oped free-con veet ion flow in a channel are (eq. (55) 
is assumed to apply with ( o— 1) 





where II is the rigid side of the [first of equations 
(61a) or (61b). 

This then is the solution of a free-conveetion 
flow wherein large temperature variations across 
the channel air admitted. The effect of longitu- 
dinal wall temperature variations would presum- 
ably be allowed for by considering that these were 
local profiles, by a sort of strip theory. Results 
of the kind in reference 4 would apply to ease III, 
discussed later on, wherein small temperature vari- 
ations are assumed throughout. 

CASE II— Tu VARIABLE, FORCED CONVECTION 

This is a flow in which the body force is con- 
sidered to be negligible. In such a ease equation 
(53) no longer applies, but one assumes that /=() 
and that e is given. Without loss, take <J r (x) = 1 - 
Equation (55) is replaced by 


1 is used here and also in equation (55) be- 
cause C 2 was introduced only as an aid in deter- 
mining what matters in the equations of motion. 
With these small terms eliminated, (\ can he 
dropped. 

It should he remembered that these equations 
correspond to oases where the temperature varia- 
tions across the channel (ran be large. The sol u- 
tions are valid provided only that bu 1 2 max jjL$T is 
small. For consistency, of course, g x (x) must be 
constant. 

The system is readily solved by first integrating 


Then equations ((H) and (62) with their accom- 
panying conditions hold, but with the right sides 
of the first of equations (61a) and (61b) replaced 
by unity. The solution is that given in equations 
(62) but with / 1 = 1 . 

The analysis thus far given applies to the ease 
where large temperature differences an* allowed 
(7^1). If only small differences are permitted, 
a somewhat different formulation results. 
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CASE III— 7^-1, FREE CONVECTION 

If T {) =\, then po— Mo=l- Also, equation (53) 
still applies, and hence the driving term in equa- 
tion (49) is (compare with eq. (54)) 


Then equations (68) and (69) arc 

■«o.»+0— 7,+! !i(s) + l+! l^ X) (70) 

5 5/37 


/. 


<Jr) ■>') 1 ■( 5p| ) 


pfi 

_ o 


I ((«ns) y (04) 
- p/S/87' — 7\ + * * J (Liquid) 

Then, fora free-eonvoetioii How, it is required that 


fc(/)=-i+o(wrj 


(65) 



(71) 


where again the terms involving e/5 in equation 
(70) and (jj. in equation (71) appear only for a gas. 

Xow to determine what form the temperature 
variations can take. It was assumed at the outset 
that u 0 is independent of x. For this case it is 
readily shown that the most general forms allow- 
able for the temperature and pressure gradient are 


ITnder such a condition the driving and viscous 
forces are of the same order provided (compare 
with eq. (55)) 


pfdfi 7 — - p ujd~ (’60) 

If these last equations and equations (51) and (52) 
are used, equations (48), (49), and (50) become 


PI'zWq 


(67) 


"(>. uu 


!\ — 


tij(j') I 
8 


-( (/iW_ »' 

8pT ff-W 


(?/„(/,, j. 


i'i "<>.„) 
(OS ) 


7 V 


W- M<V f r r | /* 

btiTr k\ ' r+ 7icr 



8 






y 


(69) 


where the irrins involving e/5 do not appear if the 
fluid is a liquid. 

Tn distinction to the result for the case of largo 
temperature variations, in equation (68), 


— ■-=:= 06, 0005 2 (For water) 

fL&T 

0,4005 2 (For air) 


/ 1 — (t\x T (i 2 -f- .7 a ( y) (72) 

(f r (x) - - 1 +fi J lJ2 7 [(e/5 + tf 1 )jr + w. 2 ] (73) 

Then equations (08) and (09) become 



n 

■0, v\t • 

- T, 



(7( 

I L\ uu 

: tttA —Wii ir 

■ 'll) 

1 t a i 

r '.( i 

. /A 

(1) 


Sk i \ 


17 A 

’ per) 

(71 

here 

lie term pH* 

'■W'{ 

Wolf 1' 

} disappears 

for 


liquid. Those equations require small temperature 
variati( ns but admit substantial mass-flow rates. 

Then , in terms of u 0 and |7 \(y), equations (70a) 
and (71 a ) are a ])a i r of ordinary differential equa- 
tions aid are nonlinear only if t he frictional heat ing 
is important. Solutions in the linear case are quite 
simple ref. 4) and in the nonlinear ease can be 
found by the iterative method given earlier follow- 
ing equ it ion (37). Ostrach discusses this system 
extensively in reference 4, where, among other 
tilings, some machine solutions art 1 given. 

CASE IV To£El, FORCED CONVECTION 

For -ase IV equation (63) still applies, arid 
(filiations (70) and (71) are replaced by 


for the circumstances of equation (59). In the 
ease for air, if <5 — e, this is a very small number 
and should, indeed, be negligible. The convection 
t erm in equation (69) is of order 5 or larger, as is t he 
frictional heating term, and should be retained. 


Wo = 0x00 (74) 



( 75 ) 
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The temperature variation must have the form 


then the permissible separated forms are 


1\ — (d-i + a 2 y ) x + T 2 (y) 


n/ui (Y) — (X/L) 6 or- e ex/L 


while g x (x) is a constant. These equations are 
easily solved because the dependent, variables are 
separated. The result is, of course, the familiar 
Poiseuille flow. 


OTHER GAS FLOWS 

The flows discussed thus far have ineludc-d only 
those cases wherein the velocity is essentially 
independent of distance along the channel. If 
this restriction is lifted, the problem becomes 
vastly more complicated. Therefore, only one 
<*lass of solutions is examined here. These solu- 
tions are ones in which t lie independent variables 
are. separable. Such a flow is out of the question 
for liquids unless if can be assumed that the 
temperature is a function of y only. This is 
becauseof the form of the state equation (eq. (5b)). 
It can be shown that this limitation to '/’independ- 
ent of x leaves only the fully developed eases 
discussed previously. 

First recall that, to have a fully developed flow 
in any sense, the channel must be very long and 
the dependence on x must he much weaker than 
that on ?/. Then equa tion (1) to (5) can be approxi- 


mated as 

(pu) x +(pv)y= 0 (76) 

P (mix ~h y ) + ? J x = — p/T (yM r) y (77) 

I\ = 0 (78) 

pc v (ul x-\-vT y ) ■F/ > (Ua’H“2V) = (kT Y ) y j tP-u\ (79) 
P—pRT (Gas) (80a) 

P=P [l— T)} (Liquid) (80b) 


These equations can be derived formally in the 
same manner as they are derived for external 
boundary layers. The only difference is that the 
Reynolds number of boundary-layer analysis is 
replaced by a ratio ( Lid ) 2 , where L is a charac- 
teristic flow length and d is the channel width. 

If the fluid is a gas (eq. (80a)) and the body 
force is negligible, and the viscosity and thermal 
conductivity each vary with the temperature as 

p = aTt (81) 

k = bTt 


r/Pi( Y) = (X/L) e or 

P / Pl (Y) = (X/L)* «r- l >+ 1 or 

T/T 1 (Y) = (X/Ly e or 

I>II\ = (X/L) 9 *t+ l '+ l or 


f 8XiL 

Y8X/L 

^ bx/l 

e <2{ + l)0X/L 


V (82) 


where 0 is an arbitrary constant. 

It is interesting, though perhaps irrelevant, that 
exactly the same variations of freo-stroam velocity 
are allowed for similar solutions of the external- 
boundary-layer equations (ref. 7). 

The exponential form in equations (82) is valid 
for the complete Na vier-Stokes equations, while 
the other form depends critically on the assump- 
tion of a very long channel. 

Several other somewhat unrelated comments 
about this result are perhaps in order. First, for 
a liquid the requirement- previously stated, that 
d r Y/&X=(), leads to the condition that nothing 
varies with X. This case has already been ex- 
amined. Second, for a gas, if the viscosity and 
conductivity do not. have the same variations 
with temperature (eqs. (81)), only the trivial A"- 
independent separation results. The A-indepetid- 
ent solution corresponds in equations (82) to the 
exponential variation with 0=0 and was discussed 
starting with equations (7). Finally, if the body 
force is important ( — pf in eq. (77)), the forms 
given in equations (82) apply, but only with 0=0 
(exponential and uninteresting) or 6~-% (power 
of X). 

The forms given in equations (82) have two 
other properties of interest.. The through-flow 
Mach number, which is proportional to u/^fT f 
is independent, of X. Also, unless there is flow 
through the channel walls, all the solutions except 
0~ 0 (exponential) and 0 — Kf (power of A') 

must be flows with zero not mass flow. This is 
because the mass flow is 


JVud r= (X/Lyw j’p, (Y )u, (Y) dr 'j 
=«***■ "j d p l (Y)u l {Y)AY 


y (8:i) 


which must not vary with X unless there is flow 
through the walls. 
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If equations (82) are put into equations (76) to 
(80a), the result is 


+ ( Pi t'i )\ — 0 


Pi 


V /: ) PlU 

/0 0 \ , r<9(2f+ 1 )-rV"| r> 

(x ■'j + L l I 1 ' 

= (p i W'i . r) y PifPi, y — 6 
"rPl (j U l T 1 + V i r i ,y )+I\( 


L 


Mrt f’l.r 


(84) 


-- (^1 7\ , Y ) Y + Pl^l, yP 1 — Pl^l / 1 

where X 0 if the exponential variation is used and 
X 1 if the variation is as a power of x. Also, the 
term p x j can appear only if, as already men- 
tioned, 0 = 0 (exponential) or 0 = 1/2 (power of ./*); 
k'x and mi have the obvious definitions. 

Now a new space variable, y, is introduced from 
equation (10). Then equations (84) become 

f PP l (2('0-|-X) ^1 , rpi _ f , \ n / cr v 

\2A'/47' t (-l)K/ Ml+r ' (p,p,) ^° (8o) 

(wt) «?+ { — -/r^-} r ' - 


+ {p 2f ^ ,)±X ] />1 } T] 


4«[7’,(-l)F ^ 

7i 2 


7’ i r M 

1 1 “I, „ y> 


(80) 






+“ / /l7’,(-l)] f W4)7’!- f 7V 

46[7’ l (— l)l 2f 

- 


-ri-f (r,. w +J«u) (s?) 
Pi=p i ItT l (88) 


Observe that equations (85), (86), and (87) form 
a set of three ordinary differential equations for 
three variables (u u T u and pith). One curious 
feat ure of the system is that it is of fifth order and 
there are, in general, six boundary conditions to 


be applied (one each to u ly v lt and T x at each wall). 
Hence some restriction must be placed on the 
combinations of boundary conditions for given 
values of the parameters, particularly I\ and 0. 
Actually this is no different from what happens in 
analysis of the external laminar boundary layer. 
In that ease, however, the boundary condition on 
v (or z’j at the outside of the boundary layer is not 
satisfied, nor is there any particular reason for it 
to be. For channel flows such an omission is 
probably not allowable. In general, this means 
simply that there must be a flow through at least 
one wall, and this flow cannot be prescribed if 
similarity is to be maintained. 

Tf, on the other hand, one examines the case 
where 2f0 + X=O (eq. (85)), it is permissible to set 
■/’! — (). Then the system is of fourth order with 
four boundary conditions. The condition 2f04 
X = () corresponds exactly to the two cases men- 
tioned in connection with equation (83). These 
are tin only ones in which i\ can vanish at the 
walls. Of the two cases defined by 2f0+X=O, 
namely, X = 0 , 0 = 0 (exponential x- variation) and 
X=l, j*0= — 1/2 (^-variation), the first has been 
solved earlier. In the second case, for linear 
viscosity-temperature variation, equations (80) 
and (87) become ( — 20=f=X= -1) 

«?+ /r/’,= - R>^7u(- i (89) 

+ ( j r, u foo) 


where he boundary conditions are 


Ui(±l)=0 

Ti(— l) = Constant 

7\(-f- 1 ) — Another constant ^ 


(91) 


The sc lut ion of this system is more difficult to 
obtain than is the solution of equations (35) to 
(37), although a similar procedure can be followed. 
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If the flow is fairly slow, with the Mach number 
limited to, say, a few tenths, the profiles are, to 
good accuracy 


[7\(+l) + T 1 (-l)] +| n\(+V)~ 7U-1)] 


Ui 


r /vp A 

\oA[7’,(+l) + 7*,(-l)]/ 

V / 1 r-y . rr. (+ 1 ) - 7’, (-1)1 *(i • —n*) Y I 

1 4 + Lr.(+ij+r,(-ijJ 12 r f J 


(92) 


where use lias been made of equation (10) to 
determine B. It may be observed that the dis- 
tance L can be defined by 



which follows directly from the last of equations 
(82). If the wall temperatures are equal, the 
velocity profile is exactly the familiar Poiseuille 
one. 

A final remark: It can be seen that the system 
described by equations (90) and (91) will probably 
admit pairs of solutions just as the free-conveetion 
flow of equations (25) to (27) does. 


CONCLUDING REMARKS 

When fully developed channel flows are con- 
sidered, the cases that can be solved exactly are 
very limited. For a gas a constant pressure is 
required, or at least one which does not vary in a 
streamwise direction. This case is analogous to 
Couette flow in that no approximations need be 
made in arriving at a relatively simple mathe- 
matical problem. In the case of a liquid, one can 
solve the exact case of constant pressure gradient 
in the streamwise direction. For both the gas 
and the liquid, the wall temperatures must be 
constant. In the gas case nothing astonishing 
happens. However, in the liquid case a sur- 
prising result arises. There appear to be (except 
for certain singular cases) either two or no solu- 
tions for the flow. This result, which has been 


discussed extensively by Ostraeh, has one solution 
for which frictional heating is negligibly small. 
The second is one in which the frictional heating is 
large, and thus the temperature is raised and the 
buoyancy effect is increased. In the present 
report an approximate analytic solution of this 
problem is given. The results agree very well 
with Ostraeh’s machine calculations. Although 
an analysis is given for the case of variable vis- 
cosity, it turns out that for the cases of interest, 
wherein the fluid velocity is kept within reason, 
the temperature variations are small and there is 
no reason to consider variable viscosity or con- 
ductivity. 

These so-called (exact solutions, particularly in 
the ease of a gas, do not cover all the flows of 
interest. Hence, consideration is given to cases 
in which there are streamwise temperature and 
pressure gradients but in which the flow velocity 
is virtually independent of distance along the 
channel. For a gas the mere presence of a pres- 
sure gradient requires a temperature gradient, 
while for a liquid the presence or absence of a 
streamwise temperature variation is governed by 
the wall temperature conditions. 

In these cases one of two situations occurs. If 
the temperature variation across the channel is of 
the order of the temperature level, then in order 
that “channel flow 7 ' he maintained, the convection 
terms in the energy equation must he negligible. 
This implies that the mean flow Mach number is 
small and also that the frictional heating is negli- 
gible. For such circumstances the equations are 
separated and can readily be integrated for any 
case of interest. Only one solution exists. 

On the other hand, if the temperature variations 
are small, more complicated effects occur. This 
situation of very small temperature changes ad- 
mits very large flow velocities (see the discussion 
following eqs. (46)). Then both the frictional 
heating and the thermal convection effects can be 
significant. In such cases (when the frictional 
heating matters) two solutions can occur. These 
flows qualify as quasi-incompressible in that the 
only place where compressibility effects matter is 
in the buoyancy term in the streamwise momen- 
tum equation. 
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The forced flow perhaps deserves an added 
remark. When the temperature varies only 
slightly across the channel, the velocity profile 
must he the usual parabolic: one. When the 
variation in temperature is large, the profiles can 
still be found in closed form but are more compli- 
cated. 

If streamwise variations of velocity are allowed, 
the flow is more complex. A description is given 
of the circumstances under which the independent 
variables are separable. These forms can yield 
new results only for gas flows and show that the 
streamwise variation must be either as a power of 
x (streamwise coordinate) or exponentially with x. 
With two exceptions, only one of which admits 
./•-variations, those flows require that the body 


force be negligible. The exponential cases apply 
to the full Navier-Stokes equations, while the 
other ones require an expansion of the equations 
of motion in terms of the width-to-length ratio of 
the ehinnel. For all these cases, the streamwise 
Mach lumber is independent of x. All but one of 
t hese possible flows lead to difficulties with bound- 
ary conditions and require a flow through the 
walls. The lone exception has streamwise velocity 
proportional to For small flow Mach 

numbers the solution is similar to that for Pois- 
euille flow but allows for temperature variations 
across the channel. 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, June 5, J958 



APPENDIX A 


SYMBOLS 



parameters defined in eqs. (41) and 
(05) 

a,b 

parameters defined in eqs. (9) 

n 

parameter defined in eqs. (10) 

a 

parameter defined in eq. (C9) 


specific heats 

d 

wall spacing 

/./, 

body forces in A”- and Indirections, 
respectively, considered positive 
in the minus A r - and indirections 

<l(s) 

pressure perturbation, eq. (47) 

//,//. 

temperature functions, eqs. (16) or 
(22) and (Cl) 

K 

parameter defined in eq. (34) 

k 

thermal conductivity 

L 

characteristic length, eq. (47) 

m 

parameter defined in eq. (37) 

P 

pressure 

It 

gas constant, eq. (5a) 

*, T a 

parameters defined in eq. (21) 

T 

temperature 

y T * 

reference temperature defined in 
eq. (19) 

r,(\ 

dimensionless velocities in eqs. (33) 
and (Cl) 

U,P 

velocity components in X - and Y- 


d iree lions, respectively 


x,r 

(Cartesian coordinates, X being in 
the main flow direction 

•r,y 

X/L and Y/d, respectively 


parameters defined in eqs. (26) 

0 

volumetric expansion coefficient (see 

<‘q- (■ r >b) 

5,6 

small parameters introduced in 
eqs. (47) 

V 

dimensionless distance across chan- 
nel in eq. (10) 


parameters defined in eqs. (SI) and 
(82) and after eq. (84) 

P 

viscosity 

P 

density 

T 

Subscripts: 

dimensionless temperature differ- 
ence in eq. (33) 

m 

mean value corresponding to aver- 
age of wall temperatures 

w 0 

wall conditions at }—0 

Wi 

wall conditions at Y~d 

X, Y, X, y, I) 

partial derivative with respect to 
that variable 

0 , 1 

Superscripts: 

zero-order and first-order solution 
in eq. (47) 

(0), (1) 

first two approximations 

bars 

reference values in eqs. (5b) or (47) 
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APPENDIX B 


A GENERALIZATION OF EQUATIONS (7) 


If 11 transverse body force /i and a distribution 
of heat sources pQ, where Q is a constant, is in- 
cluded, equations (7) become 


e=0 

(pUy)y^pf 

/'» — Pfl | 

tkTy)y —pUl—pQ j 


(B I ) 


These equations can be solved in exactly the same 
manner as equations (7). The results are, for 
linear variation of viscosity and conductivity with 
temperature, h 0, and Q—C), 




u [Cn— l)**- (v+Ve-’+'Ie-*'] 

-<ij i 


T 


I I U>„ 


+ *7 


/ o 4, 

(V— e~ a y — (e 9 — e~ 9 )e~ 9V 
a 


(i / // Vi’ y f v 
b \2aj\ ) \2 

2 

I , e~ -’<q — (e 2 " I €~~ a ) 1 

>*]} 


(B2) 


where c = Bf t 12 RT Wu is defined by 

'= 2 -£Sf, aw w 

2 ] 

and, for /, = 0 and 

/'- Constant 

m c/ , 2 (»r— • ) 


+ 


2 

PQB 2 

8 Win, 


, /T w -T m \ a/gg/PV/^-K 
2 ^ 12 / 


(1-V) 


where 


(B3) 



\ / j^ « \ = | 


If V is positive, there is only one real root of this 
equation, that root being such that 


0< 




<i 


If ^ is sufficiently negative, there can be three real 
positive roots. 

Other solutions can readily be obtained for the 
case where neither^ nor Q vanishes or where other 
distrib ltions of heat sources occur. 
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APPENDIX C 

EXACT SOLUTION FOR LIQUID 


The problem at hand is to solve equations (27) 
and (17l>), subject to equations (8). Define 


kTi+2sII(T) / T+T a V'\ 
1 r/’M- 2sH ( T*) V T* + TJ 


I 2s _ U\-2s/k 

1 ' ~ u \ wr+Mi*~r*~+T tt , 


In terms of these variables the problem is 


(Cl) 


I'u 


[B*ptifsy!2s/icl( 1 \ 

L 4(7’*+ 7„) J\ v //,/ 


//.. 


-I ’l. 


U i(± 1) =0 

//.(± 1) = [(r„ 0 + 7V)/(T*+ 7’«)] 2 


( 02 ) 

(03) 

(04) 


where, to keep the problem from getting out of 
hand, equal wall temperatures are assumed. For 
convenience, define 


Ai = [(T W0 +Ta)l(T*+T a )Y 

4 = B ! #v2# 

" 2 4(7'*+ ? a) 


(C5) 


It should be remembered that the parameter B, 
and hence is as yet undetermined (recall e<). 
( 10 )). 

The solution is found in exactly the same man- 
ner as the solution of equations (;5f>), (MO), and 
(87). Thus assume 

r (0) .1(1-+') 


Then equations (08) and ((’4) yield 

IIi=Ai +/1 2 (1 — ij 4 )/d (CO) 


If this is put into equation (02), that expression 
becomes 

vVA \ 

((7j 

\ A ' 2 v ) 


i 


If this is integrated and the appropriate boundary 
conditions are satisfied, the result- is 

...r> +f ^ \ (c 

JoVC* — y* J v \ ( 4 — y A J 

C^l+ZAJA 2 


where 


(C8) 

(CO) 


If equations (C8) and (CG) are put back into 
equation (C3) and the result is integrated across 
the channel, the result is 

+ 1 | :t a UnO' e J 


+ 


Before using this equation to determine A, some- 
thing has to be done about B, which is as yet 
undetermined. If equations (21), (Cl), and (Ob) 
are used, the second of (-((nations (10) yields 


iiA(T*+T a ) d v 


JtsB 


r dT? 
Joy+W 


(Oil) 
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If this is put into the first of equations (Co), 
eliniinating B , and the result is put into equation 
((MO), the defining equation for A is, finally, 

A'-tAyl , 3 ( 1( JjzX 2^ v - + 9 ( V I dv ' Yd 7, 

Jo V t: 4 -, 4 Jo\Jov6 ,4 -t,v 

288 S /k[Y- dv V 
- JAv^-W =0 (C12) 

[p/*W* + 7V)]' J 

This can lie solved to very good accuracy by 

r <>v v 


, ! .. . 2 V 72s/ 1 ,) 

I . 3 V 3 f ( i — r)dt? yo \'( —v / 

2 Jo vr 4 -v rf/W 


+ 7«) 


(Cl a) 


(C8) aid (Cl a) ran be evaluated as follows 
(ref. 8) ; 

p . ,iw-2,r ?(*,,-• i 

Jo VC'-V* c y 2 1_ v \r+6 K | 

— /Yir/2 , l/V-) — 1\ cos' 1 ) ?/C, 1 /y’2) J 
I 1 A ,I L =1 [siii-Hl/C^-siir 1 ! ^/C*)] 

J, A r 4 - , 4 2 1 ' /J 

-o 7 [sin - V2/( I 1 (■ , l/v'2']- | 

(Cl 4) 


The second root of equation (Cl 2) is never a 
significant one. This expiation can readily be 
solved for A in terms of ( r (which is itself a 
function of A] see eq. (C9)) and the parameter 

[rff<P(T*+T a ) ] 2 = TiT 

i (^A i — 1 j 

where K is defined in equation (84). As in the 
case of constant properties, two possible values 
for A are again found. 

The various integrals appearing in equations 


when* F(<pj] j\ ! 2) and K{ip, l/\2) are the respective 
elliptic integrals of the first and second kind of 
amplitude 1 and modulus l/y'2. Two of these 
integrals can be approximated as follows: 


d?7 _ 1 


(v+t/10 C 4 ) 


which is correct to 1 percent if (01.3, and 


f 1 3 (1— n)di» 2 

Jo N L He* 



3 

706* 



which is correct to 1 percent- if C> 1.0. 
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